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Abstract 

We compute the ratio of the coefficient of shear viscosity to entropy density at finite 
coupling and at zero chemical potential using holographic duality, up to ten derivative terms 
in the low energy effective 5-dimensional action, of a specific kind, which may or may not be 
connected to the supersymmetric completion of Type IIB theory. The result suggests that 
this ratio can be positive only for the 8th derivative term even with the form of that term 
in the action as C tj k iC kl mn C mn rs C rs ij, where C is the Weyl tensor. 
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1 Introduction 



It has become very interesting to understand the properties of the newly discovered state 
of matter, which is characterized by having very high energy density, strong collective flow 
and early thermalization time [1], at RHIC and may be called as the "first signature" for the 
formation of quark-gluon plasma. However, the strong criteria to see such a phase would be 
the quark deconfinement and chiral symmetry restoration [1]. Moreover, since the phase is 
seen at the strongly coupled limit means the quantitative study of this phase is going to be 
very cumbersome. However, we can understand the properties of some other theoretically 
constructed phases of similar kind, which may share some of its properties in common with 
sQGP. Let us assume that is the case and denote this theoretically constructed phase as 
hQGP. This means some universal properties could exist which are in common to both 
sQGP as well as to hQGP, i.e. those properties that remain the same irrespective of the 
model understudy. Now, we can understand quantitatively the properties of these phase 
using the holographic prescription [2]. 

One such quantity is the ratio of shear viscosity to entropy density, rj/s. This ratio has 
been computed earlier for M = 4 system at finite temperature [3] and is found to have the 
value j-. It has been conjectured in [4] that for any (relativistic) system that admits a gravity 
dual should have a minimum value of j- at zero chemical potential. This conjecture was 
based mainly on the calculations done (during that time) in a very specific limit that is in the 
large rank of the gauge group and for large value of the 't Hooft coupling (N — > oo, A — > oo) 
and also in the absence of enough experimental data. However, recently it has been shown 
explicitly that for some theoretically constructed models this ratio can even go below j- and 
this happens at finite 't Hooft coupling limit [5] and [6]. In [7], a theoretical argument was 
provided with the help of a Gedanken type experiments to falsify the KSS conjecture. See 
[8], for nice reviews on the subject. 

In this note, we shall present results based on the calculations done using holographic 
duality. We evaluate the ratio of rj/s up to tenth orders in derivative expansion using 
the elegant approach proposed by Iqbal and Liu [9] and is beautifully presented to study 
some more examples in [10]. In this approach [9], one important property is: the transport 
properties at the boundary is fully captured at the horizon but only in the zero frequency 
limit. This approach has the drawback, we can evaluate the transport quantities only in the 
scalar channel. 

In our computation, the only degrees of freedom we have considered is the metric. It is 
assumed not to couple with any other degrees of freedom or with itself in a manner other 
than those that are governed by Weyl tensor, which means no derivative occurs to tensors 
like Ricci, Riemann tensors and Ricci scalar. This suggests the action will be constructed 
using the Weyl tensor and we choose the contraction of indices to be of a very specific kind. 

The result of the calculation is that at the fourth power to Weyl tensor, of the type 
C lj kiC kl mn C mn r S C rs ij, we can see there is an enhancement to both the shear viscosity and the 
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entropy density and making the ratio of rj/s to increase above 1/4-7T, provided we take the 
corresponding coupling A4 to be a positive number. 

The paper is organized as follows: in section 2, we shall write down a low energy effective 
action in 5-dimension, whose rj/s we want to calculate and review briefly the computation of 
i] following [9] and entropy density using the Wald's prescription [11]. Then we will present 
the result of the computation in section 3. We have relegated the explicit computational 
details of both 77 and entropy density for two examples, the Gauss-Bonnet term and square 
of Weyl tensor in section 6 of Appendix A and the relation between A2 with the central 
charges in section 7 of Appendix B. 

The computation to rj/s for different examples both with and without the higher deriva- 
tive corrections have been studied in [10], [12], [13], [15], [16], [17], [18], [19], [20], [21], [22] 
and 



2 The action and the prescription 



Let us assume the low energy effective action that we are interested in has the following 
structure 
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where C abcd is the Weyl tensor. In order to evaluate the quantities of interest, we shall follow 
[9], which says that the shear viscosity can be evaluated by computing the following quantity 

,. U(r,k a ) 

7] = hm — - — — (2) 

ka->-o t(j(p{r, k a ) 

at the boundary, where n is the momentum associated to the field and is related to the 
metric fluctuation, in particular, to the fluctuation in scalar channel. Let us denote the 
metric fluctuation 

h\ = / [dk]Mr)e~ Mz , (3) 



where the graviton is moving along z direction and we are using a short hand notation to 
write the appropriate measure factor for momentum integrals and factors of 27i in [dk\. 

Considering this kind of fluctuation, one can show that the field <f> decouples from the rest 
of the metric fluctuation [24] and its equation of motion can be derived from the following 
effective action 

f dr[dk] ( A{r)<t>l<t>- k + B{r)^_ k + C(r)^- k + D{r)4> k 4>- k + 
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where K, is the appropriately generalized Gibbons-Hawking boundary term, which is written 
in [12], [10]. One of the important point is that upon inclusion of higher derivative terms to 
action, it generically gives us a non-trivial form of the functions E(r) and F(r). For example, 
in the 2-derivative action, like the Einstein-Hilbert type, these two functions vanish and are 
non-trivial for other higher derivative terms like Gauss-Bonnet etc. This means these two 
functions E(r) and F(r) are already of order Aj. As we know via the Kubo formula, the 
coefficient of shear viscosity is related to the two point correlation function involving the 
energy momentum tensor and our interest is to compute it to leading order in the parameters 
Aj, means we shall compute the functions A, B, ■ ■ • , F, E etc., to linear in Aj and not like 
the product of Aj's or powers of more than unity. 

By considering the radial coordinate as time, as is considered in the paper of Liu et al., 
we can compute the momentum and it comes out as 



8ttG 



B-A 
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(5) 



In order to compute the momentum in zero frequency limit and then evaluating the 
ratio eq(2) at the boundary requires us to know the derivatives of the field 0. In the zero 
frequency limit the momentum is constant [9], which means it has the same value either on 
the boundary or at the horizon. On evaluating it at the horizon with the in falling boundary 
condition for the field requires us to use 
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which is argued to work only at the horizon in [9] and we have assumed the five-dimensional 
metric has the following form, where there is a rotational invariance in the spatial directions 



ds 2 



-a(r)dt 2 + b(r)dr 2 + c{r)[dx 2 + dy 2 + dz 2 



(7) 



Using these ingredients that eq(2) gives us, after dropping the index k, of course the 
momentum dependencies are there in the field, 
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In order to evaluate the ratio i]/s, we need to know the expression to entropy density s. 
The entropy of the gravitational system can be calculated using Wald's formula [11] 
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(10) 
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where L is the Lagrangian that follows from action eq(l) and the binormal quantity e a b is 
normalized to obey e a be ab = —2. One particular choice to construct such an object is to 
introduce two null vectors £ a and 5b with the restriction £.5 = 1 on the horizon [25] 



e a 6 — (, a Sb — 8 a £,b, 
with the choice to our metric as written in eq(7), it gives 



6 = -a, ^ = 1, St = l, 5 l = 
and the rest of the components of £ a and 5b vanishes. 
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3 Computation 

Before doing the computation of r\ for non-zero but small positive \, let us write down the 
solution to metric for the case Aj = 0. This solution can be read out very easily from the 
paper [5] and is given by 
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Let us set N = 1 or can be re-absorbed in defining the temporal coordinate. Now, to do 
the computation for both r\ and S, we do not need to know the Aj corrected solution 1 . This 
is due to the reason that we are only interested to calculate these quantities at the leading 
order in Aj only and not in their product. This also means that we can calculate each of the 
higher derivative terms without the need to worry about the other terms in the action. It 
means the action that we need to consider are the Einstein-Hilbert term, the cosmological 
constant term and any one of the higher derivative term to find the result at that order in 
derivative expansion. For example, the action that we shall start out with contains only the 
A2 term 



1 
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(14) 



Going through the above mentioned procedure as outlined in section 2, this gives us the 
result that there is an enhancement to both the coefficient of shear viscosity and to the 
entropy density ( assuming Aj is positive, which we do not claim it should always be. ) 



rj 



16itGL 3 



1 + 



4A2 
L 2 



s = 



_S_ = f rj \ 
~~ V* \AGL 3 J 



(15) 



1 This we demonstrate by going through two different examples in the appendix. The main point is that 
the the uncorrected piece to both i] and entropy density depends only on the spatial part of the metric 
component, which do not receive any correction at the linear order in Aj. 
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giving the ratio 

1 = 1. 

S 47T 

We do not write down explicitly the form of k 2 and k 4 as these are very big expressions 
and similarly for other cases too. 

At this order, if we would have considered the Gauss-Bonnet term instead of the square 
of the Weyl tensor term then the result on both the coefficient of shear viscosity and the 
entropy density would have been different. In this case there occurs a suppression to rj and 
unchanged form of entropy density [5] 



V = 



but the ratio remains the same 
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as that considered by including the square of Weyl term. 
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Let us move to the next order in derivative expansion and consider the following action 



1 



16nG 



I 



-9 



(19) 



The computation results in the following form to the shear viscosity and the entropy 
density 
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and the ratio 
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If we would have taken another structure to the action at this order instead of considering 
the cubic power to Weyl tensor, like the one that follows from the third order term to Lovelock 
gravity 
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24:R ljkl R k ij n R n i + 3RR %: i kl R kh j + 2AR tjkl R ki Rij + 16R tj Rj k R k i — 

12RR ij Rij + R 3 ^j , (22) 

then the shear viscosity is same as that of computing from a 2-derivative action i.e the 
Einstein-Hilbert action. This happens because the momentum associated to the metric fluc- 
tuation vanishes, when evaluated at the horizon. 
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Let us proceed further and do the calculation by including yet another term to action, 
which is higher than the previous one that is the fourth power of Weyl tensor. The explicit 
form the action is 



S,= 



16nG 
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hi 



Cmn nrs 



The computation to shear viscosity and entropy density at this order results 
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So, we see that for the fourth power to Weyl tensor both the coefficient of shear viscosity 
and the entropy density increases only for positive A 4 . Thus making the ratio to increase 
and respects the KSS bound. 

In literature the computation at this order is done by including the other kind of contrac- 
tion of indices [12] and the result is that ratio rj/s respects the KSS bound. Here we see the 
result of respecting the KSS bound follows just by considering only one kind of contraction 
to indices. 

Let us proceed further and include the next order term that is a fifth power of Weyl 
tensor, with the structure to action 
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The computation results at this order to shear viscosity and entropy density as 
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(26) 



(27) 



(28) 



From the result it just follows that even though there is an enhancement to entropy 
density but the suppression to shear viscosity makes the ratio rj/s to go below the KSS 
bound, with the assumption that A5 is positive. 

Now we can include all the independent contributions and write the ratio to rj/s as 



77 
s 



1 

An 



A 2 A 3 _ A 4 35840 A 5 
1 - 8- - 384^ + 986--^^ 



(29) 
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which we can always do because we are interested to compute the ratio only to leading order 
in Aj and not their product. 

We can compare our result with the one that appeared recently [13] and the result up to 
the linear order in A 3 matches and the coefficient of A 4 and A 5 are new. In the paper [12], 
the authors had done the calculation by considering a C 4 term in the action but the way the 
indices are contracted and the number of such terms are different than we have considered 
in the present paper. The reason of considering such a single term at the C 4 order is just to 
show that even with one term we can get a positive contribution to the ratio i]/s. 

From the result of the computation eq(29), it just follows that to leading in Aj, we can 
rewrite it as 



71,3 = I 



£ A 



(30) 



where t|j is a constant, which can take both positive and negative values but the magnitude 
of it increases with the number of derivatives that we are considering in the low energy 
effective action. From dimensional analysis, it just follows that A« ~ a'**" 1 -*, where a' = I 2 , 
the square of string length, which by AdS/CFT duality means Aj ~ —{hry, where A is the 't 

A 2 

Hooft coupling. 

As a check to \ 2 and we find that our computation matches with the ones reported in 
[13]. 



4 Conclusion and discussion 

Using the approach of [9], we have calculated one of the transport quantity that is rj/s 
and have found some agreement with the results mentioned in [13] and predicted a couple 
more. In particular, we have shown that the KSS bound is respected by considering one 
particular type of term at the fourth power to Weyl tensor, provided we take the coupling 
A 4 as a positive quantity. To the fifth power of Weyl tensor, the ratio rj/s do not respect the 
KSS bound, once again for positive coupling A 5 . It is expected on general grounds that the 
couplings Aj in the gravity side are related to the central charges a and c on the dual field 
theory but the precise relation is not known for all Aj, at present. However, A2 is related to 
a and c via eq(50) 

— (31) 
L 2 ~ 8 a ' [6L) 

which are further related to two parameters t 2 and £4 defined in [26] with some constraints 

on the parameters t 2 and £4 coming from the argument that there should be only positive 

energy that is deposited on the calorimeter "experiment". There exists various restrictions 

on a/c, depending on the amount of supersymmetry preserved by the system [26]. For a 

system that does not preserve any amount of supersymmetry, one gets the restriction on A 2 

as < £f < \, for M = 1, it is < |f < § and for M = 2, it is < ^ < \. 
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Certainly, it is interesting to use the criteria of micro-causality violation [14] to find, 
if there exists any other constraint on the Aj's. But we do not have the explicit gravity 
solutions for all these cases, which deserve further investigations and possible violation to 
KSS bound[6]. 

One of the important question that arises in the study is the convergence of eq(30) i.e 
how small are the Aj's such that the sum in eq(30) converges ? Is the convergence going to 
be a criteria to fix the nature of the action other than supersymmetry? 

The choice of taking, the kind of action as written in eq(l) may not have been inspired 
from string theoretic suggestions, especially the 8th derivative term. If we take the string 
theoretic suggestions at this order then we already know the answer to rj/s [12] and the aim 
here is not to repeat the calculation. Rather to show that with a specific kind of contraction 
to indices at C 4 order in derivative expansion, we can still make the KSS bound to hold 
provided the coupling A 4 is positive. Of course, a priori there is not any reason to accept 
this kind of low energy effective action. It could be, one may turn around and ask the 
question: Is it possible to construct any other type of action at this order in the derivative 
expansion, such that it respect KSS bound ? And here is an example. Also, we are doing 
phenomenology and trying to find the consequences, if the low energy effective action has a 
form the kind eq(l). 
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6 Appendix A 

In this appendix, we shall demonstrate why we do not need to know the exact form of the 
Xi corrected solution, by going through two different examples, as stated in the main text. 
First we shall show for the case, when the gravitational action has two parts: one is Einstein- 
Hilbert action and the other is of the Gauss-Bonnet kind. Then we shall consider adding the 
Weyl-squared term. The result is that as long as the spatial part of the metric components 
do not receive any corrections to linear in A«, the computation of i]/s do not require the A« 
corrected solution at this order. 
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6.1 Einstein-Hilbert and Gauss-Bonnet term 

In this case the action is described by 



S 



d 5 x^/—g 



R +J^ + M(R m R ijkl - 4:RijR ij + R 2 ) 
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and it admits an exact solution of the following form 



ds 2 = -r 2 a 2 (r)dt 2 + 



dr 2 



r 2 j3 2 {r) 



+ 1 2 {r){dx 2 + dy 2 + dz 2 ). 



(32) 



(33) 



We have checked explicitly that this form of the metric solves the equations of motion and 
the form of the metric components are [5] 



a 2 (r) = (3 2 (r) = — 



I x\ i 1 - 4 



, 7 2 (r) = r 2 , 



(34) 



where we have set the size of AdS radius to unity, for convenience. As reviewed in section 2, 
we can calculate the coefficient of shear- viscosity using eq(8). After completing the necessary 
computations to find the expressions for A, B, E and F, we ended up 



A 



16ttG 

B 
16%G 

F 

16nG 
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[3a 2 b 2 c 2 + \ 2 (bc 2 a' 2 + 3ac 2 a'b' - llabca'c' + 3a 2 cb'c' - 5a 2 bc' 2 



■8nGab 2 c 



0, 



(35) 



where a, b and c are the metric components that appear in the notation of eq(7). The 
explicit form of 



«2 



.3/2 



16nG 32nG 



-A 2 
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Ki = 0, 



with 



r 

4A^ 



1- \/l-8A 2 (l-4 



c = r 



(36) 



(37) 



Computing 77 using eq(8) and keeping terms to linear in A2, gives the result as written 
in eq(17). The main point here is that the term k 2 and Ki upon evaluating at the horizon, 
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have the same structure whether we use the A2 corrected solution or not. This essentially 
says that the extra terms that contribute to /t 2 (r) vanish upon evaluating at the horizon. 

The entropy density calculated using the Wald's formula give the same answer whether 
we use the A 2 corrected solution eq(33) or the uncorrected solution eq(13). This is very 
easy to convince oneself. The reason behind this is that there appears two terms after the 
differentiation in eq(10), one from differentiating the Einstein-Hilbert term and the other is 
from differentiating the Gauss-Bonnet term. The second kind of term is already linear in A2, 
and our interest is to calculate entropy to linear order, which means the uncorrected solution 
eq(13) is good enough for this term. Whereas the first term depends on the spatial part of 
the metric components 7 2 (V), which do not receive any corrections in A2, the computation 
of entropy density too does not require the A 2 corrected solution. 

6.2 Einstein-Hilbert and Weyl-squared term 

In this case the action is described by 



16nG 

The equations of motion that results 



d b x\J—g 



R+ + ^2 CijkiC i:ikl 



(38) 



1 A 2 
Rij - -9ijR - ~^9ijW - Qgij + A 2 



R 4 



ARipRf - - \VjViR - \g l0 V 2 R + ^ V 2 % 

where W = Cij k iC ljkl is the Weyl-squared term. 

This equation of motion admits the following solution 



0, 



(39) 



ds 2 = -r 2 a 2 (r)dt 2 + 



dr 2 



r 2 /3 2 (r) 



+ 7 2 (r)(rfx 2 + dy 2 + dz 2 ). 



(40) 



Generically, it is very difficult to find the exact solution, however for our purpose the 
solution to linear in A 2 is good enough and is 



a = a = 



1-A 2 (^ 



0(A 2 ) 2 , 7 2 (r)=r 2 + 0(A 2 ) 2 . 



(41) 



It is interesting to note that 7 do not receive any correction to linear in A 2 . Once, again 
going through the procedure as presented in section 2, gives 



A = 



-[96aW + A 2 (-16aWa /2 - 16a 3 bc 4 a'b' 



48aW 
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lQa 3 b 2 c 3 a'c' + 16a%c 3 b'c' + 32aWc' 2 + 32aWa" - 32a 4 6 2 c 3 c" 



fi = ^^[ 7 2aW + A 2 (12a 2 6 2 cV 2 + 12a 3 6cV6 / + 16aW- 
48a 4 6 4 c 4L V 

52a WaV - Ua 4 bc 3 b'c' + 56a 4 o 2 cV 2 + 8aWa" - 8a 4 &Vc")], 
F = A 2 [-32a 3 o 2 cV - 64a 4 6c 4 6 / + 96a VcV], 

* - 4A 2 ^. (42) 
Substituting all these terms into eq(8), gives 

r 3/2 

r? = + A 2 ^r^[bca' 2 + aca'b' - aba'd - a 2 b'c' - 2abca" + 2a 2 bc% (43) 

167rG 327rGa 2 cr 

which need to be evaluated at the horizon, r = r . The second piece of this equation is 
already linear in A 2 . From the first term, there will not be any correction at the linear 
order as c does not receive any correction at this order. The argument given in the previous 
subsection goes through for the calculation of entropy density and it gives the result as 
written in eq(15). 



7 Appendix B: a and c central charges 

Following [27], we now write down the expressions to central charges a and c. For the gravity 
action 

S = -L f d 5 x[R - 2 A + aR 2 + (3R MN R MN + iRmnklR MNK % (44) 
the central charges are 



167T 2 

a 
16^2 



L 3 r 1 1 / _ 5a _ p 7 
16 + Z 2 V ~ T ~ 2 + 4 



/t 2 
L 3 



/"be; 



1 1 / 5a /3 7 
16 + Z 2 ^ ~ Y ~ 2 ~ 4 



(45) 



where L is the size of the AdS radius. For a choice like that of the Gauss-Bonnet combination, 
namely, a = t 2 , = — 4t 2 , 7 = £2, gives 



16tt 2 16«? 



1-4*" 
L 2 



The ratio 



a 
c 



Y2% 
-4ft 



16tt 2 16«§ 
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1 - 12^ 
L 2 



3- 



4^r 

L 2 



(46) 
(47) 
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This is the exact result, however, if we want to make contact with eq(5.1) of [5], then we 
need to identify |^ = \gb/2 and to leading order in \gb, 



16tt 2 

a 
16tt2 
a 



1-4- 



A 



GB 



1 - 12 



2 

Agb 



16/t 2 



5 

L 3 



1 -4A 
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3 ^/l-4X G B 



c VI -4A GB ' 

For a choice like that of Weyl-squared combination, a = t -^ ) (3 — — |t 2 , 1 = h, gives 

1 + 8 z? 



(48 



16tt 2 16«1 
It just follows trivially that 



16tt 2 16k?' 



(49) 



t 2 1 c — a 
I 2 ~ 8 a 



(50) 
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